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Two closely related models of oriented self-avoiding walks 
(OSAWs) on a square lattice are studied. We use the pruned- 
enriched Rosenbluth method to determine numerically the 
phase diagram. Both models have three phases: a tight-spiral 
phase in which the binding of parallel steps dominates, a col- 
lapsed phase when the binding of anti-parallel steps domi- 
nates, and a free (open coil) phase. We show that the system 
features a first-order phase transition from the free phase to 
the tight-spiral phase, while both other transitions are con- 
tinuous. The location of the phases is determined accurately. 
We also study turning numbers and gamma exponents in var- 
ious regions of the phase diagram. 



I. INTRODUCTION 

Many aspects of the behavior of polymers can be de- 
scribed by self-avoiding walks on a lattice. To incorporate 
the interactions of the polymer with itself, a binding en- 
ergy e can be assigned if either two nearest-neighbour 
lattice sites are visited by the polymer (point-contact 
model) , or if two steps of the SAW are located on opposite 
sides of a plaquette of the lattice (step-contact model). 
Both models show the same qualitative behaviour. 

To describe the coil-globule ("theta") transition, it is 
sufficient to assume these interactions to be isotropic, 
but some polymers have interactions that depend on 
the spatial orientation of the polymer, for instance A- 
B polyester. Such polymers are conveniently modeled by 
oriented self-avoiding walks (OSAW) with short-ranged 
interactionpibytween steps depending on their relative 
orientationcTlD. This orientation-dependent strength of 
the interaction can be incorporated in the model by dis- 
tinguishing parallel and anti-parallel step-contacts, and 
assigning an additional binding energy ep only between 
parallel step-contacts; with point-contact energies, incor- 
porating orientation-dependence in an elegant manner is 
complicated by the ends of the polymer. 

Most research on SAWs has been based on the point- 
contact model, since it is numerically better behaved 
than the step-contact model. Most research on OSAWs 
however has been based on the step-contact model, in 
which the inclusion of the additional binding of paral- 
lel step-contacts is more natural to the model. We will 
study both models inlhis manuscript. 

Bennett- Wood et aa enumerated all configurations up 
to SAWs with a length of n = 29 and ordered them ac- 
cording to their number of parallel and anti-parallel step- 
contacts. These results showed the existence of three 
phases: a free SAW phase, a normal collapsed phase and 



a compact spiral phase. The transition from the free to 
the spiral phase was conjectured to be of first order. 

For the case e = (only binding energies between par- 
allel step-contacts), Barkema and FlesiaQ extended the 
exact enumeration of the OSAWs to length n = 34. In 
the same paper, the energy of the ground state (for pos- 
itive Cp) and its degeneracy as a function of length was 
given, and an approximation to the number of configu- 
rations C„(mp) of polymers of length n with rup parallel 
step-contacts was proposed: 



C„(l) =p„C„(0) 
C„(m) w C„(l) • exp(-(7„(m 



1)) (m>l). 
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where Pn and q„ are ri— dependent parameters. The 
partition function for e = can be constructed, and 
from this it was concluded that the transition from 
the free phase to the spiral phase at e = occurs at 
Ep = login) = 0.9i201 (where ^ = 2.638 is the growth 
constant for SAWda) and is a first order transition. We 
use units such that ksT = 1. _ 

For the step-contact model, Prellberg and DrosselH ar- 
gued that the transition from the collapsed phase to the 
spiral phase is continuous. They also argued that the 
location of the theta-transition (from the free to the col- 
lapsed phase) is independent of Cp. 

Trovato and SenoQ performed transfer matrix calcula- 
tions on the point-contact model. They also made some 
calculations for the step-contact model, but with much 
less conclusive results. They found that the transition 
from either the free or the collapsed to the spiral phase 
is probably of first order. 

In this paper, we emclpjj the pruned-enriched Rosen- 
bluth method (PERM)El'El to study the phase diagram 
of both models for two-dimensional OSAWs. The only 
deviation from the algorithm as described in the above 
references is that new steps were biased both towards 
large numbers of contacts and large absolute values of the 
turning number, with different biases in different parts of 
the phase diagram. As usual in PERM, this is corrected 
for by reweighting. 

The manuscript is organized as follows. In section || 
we study the phase transition from the free phase towards 
the spiral phase in the step-contact model. To do this, we 
determine numerically the partition function along three 
lines e =constant, for all e^. 

In section [II, we use the same technique as in section 
H to study the transition from the collapsed to the spiral 
phase, but for different values of e which are above the 
collapse energy eg. 
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The next section, section IV, studies the transitions 
for e close to the coUapse point eg. To do this, we use 
the point-contact model, since for this model finite-size 
corrections at e = eg are smaller, and simulations using 
PERM have smaller statistical fluctuations. 

In section 0, we present our numerically determined 
phase diagram for both models, discuss the nature and 
location of all phase transitions, and summarize the re- 
sults. 



II. TRANSITION FROM THE FREE PHASE TO 
THE SPIRAL PHASE 



With PERM, we measured C„(mp), the contribution 
to the partition function for walks of n steps from con- 
figurations with TTip parallel contacts. First, we did this 
at e = 0, so that we can compare our results with pre- 
vious work. In figure [I] we have plotted log(C„(mp)) as 
a function of nip for various chain lengths n. To a good 
approximation, the curves are straight lines, in agree- 
ment with the guess eq. (0). The inset shows the de- 
viation from the straight line for n = 256, by plotting 
log(C„(mp)) -I- 1.229 mp. The figure combines data ob- 
tained at several values of tp close to the transition. 




FIG. 1. Logarithm of the contribution to the partition 
function CnirUp) of configurations with rup parallel contacts, 
in the absence of anti-parallel interactions (i.e., e = 0). The 
lines correspond to chain lengths of n = 32, 64, . . . , 256 steps. 
The inset shows the deviation from a straight line, by plotting 
log(C256(T7^p)) + 1.229 JTip as a function of mp. The presence 
of two peaks indicates a first-order phase transition. 



"bump" at some number ml of parallel contacts. Looking 
at increasing lengths n of the OSAWs, we see no evidence 
that asymptotically ml/n nor ml/n — s- 1, and there- 
fore conclude that this bump will persist in the thermo- 
dynamic limit (n — *■ oo). To extract the transition tem- 
perature ep,crit('T-)j we have determined for which the 
top of the bump and the walks with no parallel contacts 
contribute equally to the partition function. The results 
are ep,crit(n)=2.11(3), 1.6(1), 1.42(4), 1.368(7), 1.315(2), 
1.278(2), 1.250(2); and 1.229(2) for n = 32, 64, . . . , 256 
respectively. 

Since the number of parallel contacts in a tight spiral 
scales as n — 4-^71, we expect corrections of order n~^/^ 
to the critical temperature. Most likely, there are also 
corrections of order n~^, and possibly other corrections. 
Assuming however that corrections of order n~^/^ are the 
leading ones, we extrapolated our values for ep.crit("-) to 
the limit n oo, and obtained for e = in the thermo- 
dynamic limit Cp^crit = 0.90(5), with the error mainly due 
to the uncertainty in the finite-size correction. 

We also studied the contribution to the partition func- 
tion as a function of the turning number t: the number 
of turns that the walk has made clockwise, minus the 
number of turns anti-clockwise. Note that the turning 
number is not equal to the wiiiding number w as de- 
fined by Duplantier and SaleuiO; for large chains in the 
free and the collapsed phase however, the two quantities 
seem to be related by ((§t)^) — 2{'w^). The factor of 
two is explained by observing that the turning number 
receives contribution from both ends of the chain, while 
only one end contributes to the winding number. For the 
spiral ground state, the turning number is roughly equal 
to 2y^{n). 




Turning number 



From the partition function plot we conclude that there 
is a first-order transition from the free phase, where con- 
figurations with no parallel contacts dominate, to the spi- 
ral phase, where configurations with many parallel con- 
tacts dominate. A closer look shows that the curves in 
figj^ are not completely straight but S-shaped, with a 



FIG. 2. Probability of turning number t as a function of 
t, for chains of length n = 256 for the case e = 0. Different 
curves correspond to tp = (dotted line), tp = 1.194 (dashed 
line), and tp = 1.253 (solid line). 
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Results for n = 256 are presented in figure H, where we 
averaged the histograms for positive and negative turning 
numbers. We observe that for ep = 1.253, the histogram 
has two peaks at ±23: we are in the spiral phase. At 
Cp = 0, the turning number is between -10 and 10, with 
a maximum at zero: we are outside the spiral phase. At 
e = 1.194 we are close to the transition, and the peaks at 
t = ±21 and the peak around zero coexist. The fact that 
in the histogram for turning numbers the peaks maintain 
their location, while their relative importance changes, 
is consistent with our earlier conclusion that the phase 
transition is first-order. For a continuous transition, we 
would expect that the two peaks in the spiral phase would 
approach zero gradually. 

We repeated this procedure for e = 0.993 and for the 
case where anti-parallel contacts are strictly forbidden 
(e — > —oo) but parallel contacts have a finite binding en- 
ergy (e ± ep is finite) . Qualitatively, the behavior is the 
same as for e = 0, and we conclude that the transition is 
also first-order. The bump seems to shift to the left with 
increasing e. For e — 0.993 we obtained in the thermo- 
dynamic limit ep^crit = 0.05(5), while if anti-parallel con- 
tacts are strictly forbidden, we find ep.crit ± e = 0.75(5). 



III. TRANSITION FROM THE COLLAPSED 
PHASE TO THE SPIRAL PHASE 

To study the transition from the collapsed phase to 
the spiral phase in the step-contact model, we used the 
same technique as in section |l[ at a particular value 
for e, we calculated the contribution of configurations to 
the partition function as a function of the number of its 
parallel contacts. This allows us to determine both the 
nature and the location of the phase transition at the 
particular value of e. 

We used this method for e =1.253, 1.435 and 1.609. 
These values are well above the theta-transition, which 
is estimated to be e = 1.21 (see section IV). Results anal- 
ogous to those shown in fig.|l], but now for e = 1.609, are 
presented in figure H There are still two maxima in the 
histogram - one at irip ~ and the other at ml > - but 
the valley between them is very shallow and much more 
narrow. Comparison of different chain lengths now sug- 
gests that in the thermodynamic limit ml/n approaches 
zero, indicating that the transition has become continu- 
ous. 

We obtained estimates for the location of the phase 
transition line that are consistent with ep.crit — 0- 




100 



FIG. 3. Logarithm of the contribution to the partition 
function Cn{mp) of configurations with nip parallel contacts, 
for the case e = 1.609 and no interaction between parallel 
steps (e+Ep ~ 0). Different curves correspond to chain lengths 
n = 32, 64, . . . , 256 steps. The inset shows the deviation from 
the straight line by plotting log(C256(m-p)) + 1-91 nip as a 
function of irip. The two peaks that were present in the case 
e = have nearly disappeared and are compatible with being 
finite size efi'ects, indicating that the first-order phase transi- 
tion has changed into a continuous one. 

The histogram of the turning numbers is plotted in fig- 
ure ^ The picture is quite different from that of figure 
^: the peaks at high positive and negative turning num- 
bers do not maintain their location if the transition line 
is approached, but shift towards zero, where they merge. 
This is consistent with a continuous phase transition. 
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FIG. 4. Probability of turning number f as a function of t, 
for chains of length n — 256, in the case e — 1.609. Different 
curves correspond to ep = 1.82 (dotted line), 1.92 (dashed 
line), and = 1.96 (solid line). 
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IV. TRANSITION FROM THE FREE TO THE 
COLLAPSED PHASE 

At first, we simulated the step-contact model at various 
values of e, to obtain a precise value of eg. Requiring that 
the end-to-end distance scales as n^/^ and, the partition 
sum as iign}/"^ (see Duplantier and SaleuiO), we got eg = 
1.21(2), independent of the value of Cp, as long as ep < 0. 
Next, we simulated the ordinary (non-oriented) point- 
contact model at various values of e, to obtain a precise 
value of eg, and we got f.o .— 0.667(1), in good agreement 
with earlier estimatesoEj. 



whether this is really the asymptotic behavior. Thus, 
while parallel bonds are unimportant below the 0-point, 
they become important above it. Consistent with this, we 
found that the probability P,i(0) — Cn{0) / J2m ^n{m) of 
having no parallel bond at all decreases to zero for e > eg, 
while it converges to a finite value for e < eg (figj^). 
Based on the analogy with parts of percolation cluster 
hulls, it was conjectured by Prellberg and DrosselQ that 
P„(0) = n"^/'' exactly at the 0-point. This is consistent 
with our data, although our data show again very slow 
convergence and would suggest an exponent « 1/4 rather 
than 2/7. 
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FIG. 5. Log-log plot of the average number of parallel 
bonds versus chain length for different values of e, and for 
ep = 0. The S-point is at e = 0.667. 

This difference in eg is due to the fact that point con- 
tacts are roughly twice as frequent as bond contacts near 
the theta point, and it makes simulations in the regime 
e > eg much harder in the step-contact model than in 
the point-contact model: due to the large value of e, 
the Boltzmann weights of different configurations fluc- 
tuate strongly, which creates problems for PERM. The 
point-contact model can be simulated more efficiently by 
PERM (error bars decrease by roughly one order of mag- 
nitude for the same CPU times), and systematic errors 
due to finite-size corrections decrease, although they stay 
sizeable in both models (the same was found by Trovato 
and SenoQ for transfer matrix calculations). 

For this reason we used the point-contact model to 
study transitions for e « eg in detail. This includes the 
coil-globule transition for ep < which happens exactly 
at e = ee, as well as the region around the triple point 
(e = eg,ep = 0). 

In the same runs we also measured the average number 
of parallel contacts. Results are shown in fig.^. They 
indicate that {nip) converges for n — > op to a finite value, 
in agreement with Barkema and FlesiaO, as long ase < eg. 
This is however no longer true for e > eg. Exactly at 
the 9-j>omt, {nip) increases roughly as for large n. 
But finite size corrections are so large that it is not clear 



FIG. 6. Log-log plot of P„{0), the chance to have no par- 
allel bonds in an n-step walk. 

The fact that P„(0) does not decrease exponentially 
with n for e < eg shows that indeed the open coil/globule 
collapse occuKai-fer all ep < at the same value of e, 
namely e = eeHu^El. Parallel contacts are simply too rare 
to effect phase boundaries for ep < 0. On the other hand, 
the fact that {nip) diverges at eg implies that this is no 
longer true for ep > 0. Thus the phase boundary has a 
singularity at (e, ep) = {eg,0), which suggests that this 
point is indeed ihe triple point where all three phase 
boundaries meetQ. 

In order to verify this and to determine the orders of 
the coil-spiral and globule-spiral transitions, we measured 
also the distribution 

P„(mp) = C„(mp)/^C„(m). (2) 

Typical results for e < eg and for e > ee are shown in 
panels (a) and (b) of fig.0, respectively. While Pn{nip) 
decreases roughly exponentially with mp in both plots, 
details are rather different. In panel (a) the exponent is 
nearly independent of n, suggesting that it is non-zero 
also for n — > oo. Thus, the free - spiral transition hap- 
pens at a positive ep. In contrast, the exponent depends 
strongly on n in pancl(b), and seems to converge to zero 
for n ^ oo. This is confirmed by a more careful analysis. 
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It shows that the collapsed - spiral transition hapgens 
exactly at ep = 0, as suggested by Trovato and SenoQ. 
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FIG. 7. Distributions P„(mp) for finding nip parallel bonds 
in chains of length n, normalized to P„(0) = 1. Each panel 
contains curves for n — 163,334,664,1307,2561, and 5005. 
Panel (a) is for e = 0.531 < ee, while panel (b) is for 
e = 0.742 > eg. 

To determine the order of the transitions, we plot 
e°'^^^"^''Pn{mp) with the parameter a(e) determined such 
that both peaks in this function have the same height 
(compare the inserts in figs.|l| and ||). This is done for 
several values of e, but only for a single chain length 
{n = 2561). Again the data are normalized to P„(0) = 1. 
Results are shown in figj8[ We see that there are two 
peaks for all values of e, but that the right peak is lo- 
cated at very small values of rup in the collapsed region, 
and moves to larger values of rup only if we go with e be- 
low the 0-point. Thus we see again that the double peak 
structure is a finite-size effect in the collapsed phase, as 
we had already seen in sec. Ill, and that the collapse-to- 
spiral transition is second order. 
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FIG. 8. Log-log plot of e''^'''"'' P„(0), with a{e) such that 
both peaks have the same height. Again normalization is such 
that P„(0) = 1. 

Finally, we also measured turning numbers at and near 
e = ee. Average squared turning numbers at the triple 
point (e, ep) = (eg, 0) and at the point (eg, —oo) are shown 
in fig.0. Apart from the by now familiar large deviations 
for small n, we see clear indications for logarithmic laws 
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compatible with-the predictions C = 24/7 at (ee, 0)113 and 
6/7 at {eg, — oo)l3 which are indicated in fig.||by straight 
lines. We have not measured turning numbers with simi- 
lar precision in other phases, but the overall picture seems 
fully compatible with that of Prellberg and DrosseE. On 



the collapsed-spiral transition line (e > eg,ep ~ 
seems to increase faster than logn (see also fig.p 
our data are less precise there. 
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FIG. 9. Average squared turning number (t^) against log n, 
for (e, ep) = (eg, 0) (+), (eg, -oo) (x), and (0.742, 0) (A). The 
straight lines are the theoretical predictions for the first two 
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V. CONCLUSIONS AND DISCUSSION OF THE 
PHASE DIAGRAM 

We have found three phases for two-dimensional OS- 
AWs: two of them, the free phase and the collapsed 
phase, also exist in normal SAWs, and have a turning 
number around zero; the third phase, the spiral phase, is 
unique for OSAWs, and has a high turning number. 

we have confirmed that the location 



In section IV 



of the transition from the free to the collapsed phase 
is independent of the strength of parallel interactions. 
The critical value is estimated to be eg = 1.21(2) for 
the step-contact model, and eg = 0.667(1) for the point- 
contact model. This transition also exists for SAWs, and 
is known to be continuous. 

The transition from the collapsed to the spiral phase 
is found to be also continuous, and located at ep = 0, in 
agreement with theoretical predictionsQ. 

In section ^ we concluded that the transition from 
the free phase to the spiral phase is a first-order one. For 
the step-contact model, we have located three points (e-|- 
ep,e) on the phase transition line: (0.75, — oo); (0.90,0), 
and (1.04, 0.993). The results for the step-contact model 
are combined in figure |l^, where the phase diagram of 
the step-contact model is presented. The phase diagram 
for the point contact model is identical except for the 
detailed location of the transition lines. 
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all decreases to zero for e > eg, while it converges to a 
finite value for e < eg. Exactly at the jtheta-point, it 
was conjectured by Prellberg and Drosselfl that P„(0) = 
^-2/7^ This is consistent with our data, although slightly 
smaller exponents are not ruled out. 

At the triple point (e, ep) = (ee,0) and at the point 
{eg, — oo), the average squared turning numbers grow log- 
artihmically with ?j— with constants as predicted|-by Du- 
plantier and Saleudlj, and Prellberg and DrosselO. 
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FIG. 10. Schematic drawing of the phase diagram of the 
step-contact model. The model has three phases, a free phase, 
a collapsed phase, and a spiral phase. The transitions from 
the collapsed phase to the free or the spiral phase are con- 
tinuous, the transition from the free to the spiral phase is of 
first order. The transition line from the collapsed to the spiral 
phase is located at tp = 0. The transition from the free to 
the collapsed phase is located at e = 1.21(2). We determined 
three points on the transition line from the free to the spiral 
phase, and connected these points as a guide to the eye. 

The probability Pn(0) of having no parallel bond at 
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